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ABSTRACT 


The computer solution of matrix equations is often difficult to 
accomplish due to an ill-conditioned matrix or high noise levels. This 
paper compares two methods of solution for matrices of various degrees 
of ill-conditioning and for various noise levels in the right-hand side 
vector. One method employs the usual Gaussian elimination. The other 
solves the equation by an optimization technique and employs a function 


minimization subroutine. 
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1 . INTRODUCTION 

It is a common occurrence in engineering and science to find the 
equation involving an unknown function expressed as an equality where 
the left-hand side consists of an operator operating on the unknown 
function and the right-hand side consists of a known function. Frequently, 
if the known and unknown functions have limited domains, then the 
method of moments may be used to express both functions as vectors and 
the operator as a matrix. The problem of solving for the unknown 
function then is reduced to a linear algebraic problem. 

Frequently, the right-hand side function is a measured quantity 
and therefore can be expected to contain some amount of noise. This 
paper compares two methods for solving a matrix equation with various 
amounts of noise in the right-hand side' vector. The first method is the 
most commonly used and employs Gaussian elimination. . The goal at the 
outset of this project was to find a method which could solve a matrix 
equation of higher condition number and greater noise level than 
previously possible using a Gaussian elimination— type subroutine. 

The second method uses an optimization technique, wherein the magnitude 
of the difference between the right- and left-hand sides is minimized. 

A matrix was found whose condition number could be systematically 
varied. At the same time varying amounts of noise were introduced into 
the right-hand side vector. The above two methods were compared, and 
finally two illustrative problems were solved using the optimization 


technique. 
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THE TEST MATRIX 
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To test the two methods , a matrix whose condition number could be 
systematically varied was necessary. The condition number used in this 
section is given by 

u 

condition number = — 

X 


8 max 
^ min 


A matrix whose off-diagonal elements are all equal to one and 
whose diagonal elements are equal to one plus delta has a simple 
expression for the condition number. It is first necessary to derive 
a formula for the determinant of such a matrix. 

For an N x N matrix called A, consider 


J N ' IV 


1 + 6 1 
1 1 + 6 


1+6 


where N stands for the dimension. Expanding J N by co-factors we obtain 


J N ' (1 + 5)J K-1 - K 2 + K 3 ' K 4 ' • • “h 


where is equal to the minor which results after deleting the j-th 
column and first row. However, each is just equal to -K^ since a 

single row interchange converts to K ^ . Thus, except for sign 

changes, all the K's are equal.. Letting we obtain 



3 


- K 2 =- 


1 

1 + 6 


1 + 6 


It can be shown by induction that the formula is 


T * N JL 

J XT = o + N6 
N 


( 1 ) 


We now use this result to find the eigenvalues of the matrix A. 
The eigenvalue equation is 


Ia - A I = 0 


( 2 ) 


Using Equations (1) and (2) , we get 


N N-l 

(6- - \y + N(6 - X) = 0 


X V X 2> ' ' ' X N-1 X min. ~ 6 


X, = X = 6 + N 
N max. 


condition number = 




max. 


X 


min. 


<5 + N 


condition number - 




(3) 



We are now able to construct a matrix of desired condition number 


merely by choosing the proper <$ and N. We next test the two methods 
for solving matrix equations. 
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3. COMPARISON OF SUBROUTINES 


The first method makes use of DGELG, an IBM SSP subroutine, 
modified to include complex numbers. DGELG solves the matrix equation 
by means of Gaussian elimination. Let A be the matrix, x the unknown 
vector, and y the known vector. Then to solve 


Ax = y 


for x, one merely supplies DGELG with A and y, and DGELG returns x. 

The second method attempts to find x by minimizing the square 
magnitude of Ax - y. That is, minimize h(x) where 

h(x) = (Ax - y)^ (Ax - y) (4) 

and + signifies Hermitian adjoint. If we let the superscripts R and I 
stand for real and imaginary parts respectively, then Equation (4) 
expands into 

h = | (Ax - y) R | 2 + | (Ax - y) 1 ! 2 

= [ (Ax - y) R ] T (Ax - y) R + [ (Ax - y) 1 ] 1 (Ax - y) 1 


/a R R . -L J_ X , I / A X R i 

= (Ax - Ax - y) + I (A x + A x - y) 


.1 I 


R% i 2 


I R 


,R I 


<a r 


-A' 


lx 


A 


Thus the N-diraensional complex problem has been reduced to a 
2N-dimensional real problem. 


A 1 , x T and y T will be defined as follows 
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A’ 




\ 


x' 



y f 



Then, 


h = (A’x’ - y * ) T (A'x* - y') 
where T signifies the transpose. Expanding h, we get 

h = (x'^A 1 ^ - y'^) (A’x' - y') 

fp m rp rri 

= x’Wx' - y’Vx* - x'V y' + y’y' 

T T T T 

= x’A' A'x' - 2x' A' y’ + y* y* 

7h = 2 A ,T A'x' - 2A ' T y ' 

= 2(A'^A'x' - A'^y') 

We now have explicit formulas for the function value and its 
gradient vector. The IBM SSP subroutine DFMFP was chosen to minimize 
this function. DFMFP uses a variable metric version of a conjugate 
gradient algorithm devised by Fletcher and Powell.'*' 
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The program used consists of one main and four subroutines. The 
first subroutine, OPER, constructs the matrix A’ and the vector y' with 
varying amounts of noise. To construct a specific y', x* was chosen to be 

“0 “ 

1 

1 

x* = 0 

0 
0 
0 

_0 _ 

The problem, then, is to recover this x', given A' and y'. 

T T 

The subroutine PRELIM constructed A r A* and A T Y’ referred to as 
ATA and ATY, respectively, in PRELIM. The subroutine EUNCT supplied 
the function value h(x) and the gradient vector for a given x. This 
method requires an initial guess be given for x' . It was chosen such 
that each element of x’ was equal to one-half. 

Random noise at percentage levels of from 0.0 to 2.25 percent was 
added to the y' vector on successive iterations by OPER. The noise was 
supplied by a random number generator yielding numbers between zero and 
one in an even distribution. 

If the magnitude of the vector, equal to the difference between the 
calculated x f vector and the exact x' vector, was less than 10 percent 
of the magnitude of the exact x* vector, the solution was considered 
successful. The results of the tests are tabulated in Table 3.1. The 
condition numbers were computed from Equation (3). 

Under no-noise conditions, DGELG was able to solve the equation 
up to the point where the condition number became equal to 10"*"^. At this 
point the diagonal elements became equal to 1.0 plus 10 ^ which con- 
tains the maximum number of significant digits that the machine can hold. 



PERCENTAGE ERROR BETWEEN THE CALCULATED AND 
TECHNIQUES UNDER NOISY CONDITIONS. G II 


Condition 

Number 

0.0 

0.25 

0.50 

: 

0.75 

510 x 

O 

O 

iH 

0.0 

G 

-.3 G 

1.1 G 

.9 G 



0.0 

0 

.40 

.9 0 

.9 0 

1.7 x 

10 1 

0.0 

G 

. .8 G 

4.7 G 

4.7 G 



0.0s 

0 

...6 0 

2.0 0 

2.4 0 

6.5 x 

10 1 

0.0 

G 

4.7 G 





0.0 

0 

3.0 0 

6.7 0 


2.6 x 

10 2 

0.0 

G 






0.0 

0 




4.1 x 

10 3 

0.0 

G 






0.0 

0 




1.0 x 

10 4 

0.0 

G 






0.0 

0 




1.0 x 

10 5 

0.0 

G 






0.0 

0 




1.0 x 

10 6 

0.0 

G 






0.0 

0 




1.0 x 

10 7 

0.0 

G 
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0.0 

0 




t 

1.0 X 

io 15 

0.0 

G 






0.0 

0 





3.1 

TRUE X VECTORS FOR GAUSSIAN AND OPTIMIZATION 
lTES GAUSSIAN. 0 INDICATES OPTIMIZATION. 


Noise 


1.00 

1.25 

1.50 

1.75 

2.00 

2,25 

o o 

1 — 1 1 — 1 
♦ ♦ 
CM CM 

1.4 G 
2.0 0 

3.2 G 
3.8 0 

2.4 G 
3.3 0 

3.7 G 
2.2 0 

5.1 G 
4.8 0 

3.6 G 
6.9 0 

3.8 G 

4.9 0 

7.2 G 


6.1 0 



c 

Percent Error = where x - the calculated x 

M 

x = the true x 


Omission implies error > 10 percent 


00 



The computer cannot hold a more Ill-conditioned matrix of the form in 
question. 

For large numbers the function value and gradient had to be 
multiplied by a large number within FUNCT. If this number was too large, 
overflow resulted within DFMFP. If this number was too small, underflow 
resulted. The number which gave the best results was 10 „ The results 

of both methods are summarized in Table 3.1. 



4. „ FIRST EXAMPLE 
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In this example a mode-matching technique is employed to solve a 

bifurcated, infinite, parallel-plane,. ..waveguide problem. The waveguide 

in question is pictured in Figure (1). If a TE^ Q wave is incident from 

the left, then the resulting fields can be determined by first solving 

the wave equation for the infinite number of normal modes in each 

region and then matching at z = 0. This results in an infinite number 

of linear equations in an infinite number of unknowns. The solution 

may be approximated by truncating the coefficient matrix and both the 

known and unknown vectors. The elements of the unknown vector correspond 

to the amplitudes of the normal modes in region A. A more complete 

discussion of the problem and the exact formulas for the matrix and 

vector elements may be found in Analytical Techniques in the Theory of 
2 

Guided Waves . 

For this example, the matrix was limited to 25 by 25. The distance 

between the outermost conductors was made equal to 25 units. The 

problem was solved for a center conductor located at distances from 

10 to 18 units from the top conductor. In the matrix generating 

program, the parameter p is equal to this distance. The subroutine 

generating the truncated matrix and known vector was written by 

3 

C. A. Klein, who attempted a solution by Gaussian elimination. His 
results are tabulated in Tables 4.1 through 4.5. 

The same problem was solved using the optimization technique with 
0.0 percent and 5.0 percent noise introduced into the known vector. 

Klein also computed the condition numbers by using the new definition 
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TABLE 4.1. NORMAL MODE AMPLITUDES FOR BIFURCATED WAVEGUIDE 

p = 10 

% Noise =0.0 
Method - Gauss 


Real(x) 

-0.693152707549 
0.155337724329 
0.317913267895D-01 
-0.733016799930D-01 
0.303396945421D-01 
0.303679488886D-01 
-0.4842 7 6679557D-01 
0. 109584438701D-01 
0. 4088109081440-01 
-0. 494418659956D-01 
-0. 543902557502D-02 
0. 778883268731D-01 
-0. 776280356439D-01 
-0. 660471391558D-01 
0.319067252346 
-0.354665663973 
-3.63453198601 
-8.28472689191 
-10.9104185003 
-9.70015620899 
-6.10127601639 
-2.71727834179 
-0.825048941140 
-0.154679367040 
-0. 136187 07 0269D-01 


Imaginary (x) 

-0.720790762992 
-0.364891591341 
-0. 746784972755D-01 
0.172187192624 
-0. 712685825042D-01 
-0.713349525603D-01 
0.113757613624 
-0.257416158227D-01 
-0. 960305447226D-01 
0.116139986216 
0.127 7 63858 138D-01 
-0.182961322905 
0.182349893318 
0.155146123165 
-0.749495706601 
0.833117126476 
8.53759230698 
19.4609982387 
25.6288032169 
22.7858715639 
14.3320260716 
6.38294414712 
1.93805736755 
0.363345096218 
0. 319906301001D-01 
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TABLE 4.2. NORMAL MODE AMPLITUDES FOR BIFURCATED WAVEGUIDE 

p * 12 

% Noise - 0:.0 
Method - Gauss 


Real(x) 

-0.612978876221 
0.158265570419 
0.2971827 92634D-01 
-0.628349839686D-01 
0. 236976248805D-01 
0.214179469823D-01 
-0. 304892338095D-01 
0. 607189332784D-02 
0.19 585881 6160D-01 
-0. 200289 933433D-01 
-0. 181000234917D-02 
0. 204861632204D-01 
-0. 152907869413D-01 
-0. 899164794132D-02 
0.2 63041303 7 98D-01 
-0.13607 9851 7 94D-01 
-0. 260088498294D-01 
0. 578118299740D-01 
-0. 164515724349D-01 
-0. 387550833S89 
-0.720974673427 
-0.677714093308 
-0.372769785929 
-0.116096347718 
-0.160859761550D-01 


Imaginary (x) 

-0.790099295853 
-0.323097392011 
-0. 606695347549D-01 
0.128276917042 
-0. 483784361655D-01 
-0. 437244992318D-01 
0.622434298386D-01 
-0.123957023223D-01 
-0. 399843582757D-01 
0.408889659114D-01 
0. 369509955324D-02 
-0. "4182227 31026D-01 
0. 312159705326D-01 
0.183563487120D-01 
-O'. 5369 9587 9920D-01 
0.277805495556D-01 
0.5309677 60506D-01 
-0.118022204340 
0. 335857011358D-01 
0.791180693118 
1.47186173261 
1.38354574211 
0.761005349015 
0.237009395475 
0. 3283934 04 191D-01 
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TABLE 4.3. NORMAL MODE AMPLITUDES FOR BIFURCATED WAVEGUIDE 


p = 14 

% Noise - 0.0 
Method - Gauss 


Real(x) 

-0.562347051281 
0.158937293771 
0.283805678513D-01 
-0. 570667538155D-01 
0. 204014765122D-01 
0.173978355082D-01 
-0. 232337332473D-01 
0. 431047696895D-02 
6. 128448798674D-01 
-0.120121195993D-01 
-0. 9804 122 31164D-03 
0.986816122027D-02 
-0.64227 6774098D-02 
-0.321036596398D-02 
0. 771365634861D-02 
-0. 31235885894 9D-02 
-0. 434970895889D-02 
0. 6263982055O1D-02 
-0. 919874992913D-03 
-0. 564189590637D-02 
0. 618814322022D-02 
0.157 94 67 5377 8D-02 
-0. 190921166762D-01 
-0. 210777001712D-01 
-0.706265865991D-02 


Imaginary (x) 

-0.826901320543 
-0.300296064468 
-0. 536222344729D-01 
0.107821903696 
-0. 385465422279D-01 
-0.328714640184D-01 
0.438978070629D-01 
-0. 814421360174D-02 
-0. 242691113031D-01 
0.226956943585D-01 
0.185239050942D-02 
-0.186449001846D-01 
0 . 1213517 50105D-01 
0.606566427305D-02 
-0. 145741794716D-01 
0,59017 0713358D-02 
0. 821833850912D-02 
-0. 118351653937D-01 
0.17 38011473 7 0D-02 
0.106597960530D-01 
-0.116918755271D-01 
-0. 298424538569D-02 
0. 360726382347D-01 
0.398241989608D-01 
0.133441846777D-01 
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TABLE 4.4. NORMAL MODE AMPLITUDES FOR BIFURCATED WAVEGUIDE 
p = 16 

% Noise = 0.0 
Method - Gauss 


Real (x) 

-0.544655824876 
0.159002204125 
0.27 9127 637 618D-01 
-0.551811109372D-01 
0. 193746450077D-01 
0.162021930016D-01 
-0. 211782419486D-01 
0.383719356428D-02 
0. 111371878738D-01 
-0.101119925885D-01 
-0. 798234279886D-03 
0. 773455715437D-02 
-0.481838950870D-02 
-0.22887 04694 74D-02 
0. 517804760418D-02 
-0. 195088330013D-02 
-0. 248733502877D-02 
0. 320720168505D-02 
-0. 408219513817D-03 
-0. 206342688285D-02 
0. 171162009022D-02 
0.278188244359D-03 
-0. 129949453829D-02 
0. 61903457124 9D-03 
0. 269179439590D-03 


Imaginary (x) 

-0.838659664243 
-0.292852620963 
-0, 514101428403D-01 
0.101633389641 
-0.356845089159D-01 
-0. 298414396956D-01 
0. 390064005475D-01 
-0 . 706740009441D-02 
-0.205126380288D-01 
0. 186244181268D-01 
0. 147019975159D-02 
-0. 142456222359D-01 
0.887458135696D-02 
0. 421537029726D-02 
-0. 953700497864D-02 
0. 359316583554D-02 
0. 458121059641D-02 
-0.590707 1 6949 6D-02 
0. 751864763181D-03 
0.380045027762D-02 
-0. 315248730212D-02 
-0. 512371239947D-03 
0. 239342833995D-02 
-0. 114014706686D-02 
-0. 495778689530D-03 
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TABLE 4.5. NORMAL MODE AMPLITUDES FOR BIFURCATED WAVEGUIDE 

p - 18 

% Noise =0.0 
Method - Gauss 


Real (x) 

-0.565127177149 
0.158936666466 
0.28463127 9896D-01 
-0.574064797457D-01 
0.205920971239D-01 
0.176275737892D-01 
-0. 236441667722D-01 
0.440910063148D-02 
0.132180244197D-01 
-0. 124499163053D-01 
-0.10249 6495929D-02 
0. 104264309693D-01 
-0. 687634469470D-02 
-0.349539382475D-02 
0.858518155290D-02 
-0. 358117554900D-02 
-0. 5199354 62 265D-02 
0. 796696919662D-02 
-0. 129329530044D-02 
-0. 954607825673D-02 
0.164964373311D-01 
-0.1587 33478035D- 01 
0.101404 6069 56D-01 
-0.4 181169 945 58D-02 
0. 882768044447D-03 


Imaginary (x) 

-0.825003802202 
-0.301521151135 
-0. 539 97830124 3D-01 
0.108906700011 
-0. 390655785550D-01 
-0. 334415365493D-01 
0. 448557059949D-01 
-0. 836457141979D-02 
-0. 250761138220D-01 
0. 236189243136D-01 
0. 194447 650923D-02 
-0. 197801397122D-01 
0.130452174067D-01 
0.663116443264D-02 
-0. 1 628707 7627 2D-01 
0. 679390223769D-02 
0. 986377420543D-02 
-0. 151142576263D-01 
0. 245353005330D-02 
0.181100092811D-01 
-0. 3129564 04857D-01 
0.301135679293D-01 
-0.192376211857D-01 
0. 793216067202D-02 
-0. 167471259380D-02 
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25 

CONDITION NUMBER = 

y I A . . [ max . 

i-i 1J 

The matrix was best conditioned at p = 16. 

Under no-noise conditions and for p greater than 10, both methods 
obtained vectors whose respective elements were within 10 percent of 
each other. This was also true for p taking on values from fourteen 
through eighteen with 5 percent noise. In the remaining noisy cases, 
the optimization method was unable to get good results. The resulting 
data using the optimization method appears in Tables 4.6 through 4.15. 


I I A | max. 
i-1 J 



TABLE 4.6. NORMAL MODE AMPLITUDES FOR BIFURCATED WAVEGUIDE 


p = 10 

% Noise =0.0 
Method - Optimization 


Real (x) 

-0. 0316590234117131600000000 
0. 1309838306768769000000000 
-0. 1188739895291309000000000 
-0.9169875059763859000000000 
-1.3313092571164340000000000 
-0. 6423843595574583000000000 
0. 5334170618741029000000000 
1. 1179309577169140000000000 
0.8934457631420527000000000 
0. 4088685933280350000000000 
0. 1065223000752136000000000 
0. 0124624326021052400000000 
-0. 6616332407926009000000000 
0. 1568329392604192000000000 
0. 0310059967205517700000000 
-0. 0690563521603437000000000 
0. 2754422488657253000000000 
0. 0264791865513247100000000 
-0. 0403854591606213300000000 
0. 0086954642420436330000000 
0.0306732203033356300000000 
-0. 0348078253225158700000000 
-0. 0035584634863269560000000 
0.0467766948546173300000000 
-0. 0421097171728256600000000 


Imag inary (x) 

0.0702629983574413000000000 
-0.2908292749065581000000000 
0. 2641218745565664000000000 
2.0398984073276960000000000 
2.9698264472788650000000000 
1.4518672054891210000000000 
-1 . 1544825515623110000000000 
-2.4583865650412570000000000- 
-1. 9706664735657690000000000 
-0.9030061507109290000000000 
-0.2354325683404167000000000 
-0. 0275575878979049900000000 
-0. 7497989729616068000000000 
-0. 3475723585230524000000000 
-0 . 0687186475916295000000000 
0. 1530573218392616000000000 
-0. 0610526803673600000000000 
-0.0586955945822449000000000 
0. 0895274885290711000000000 
-0. 0192778921910766000000000 
-0.0680092264481303000000000 
0. 0771854266295259000000000 
0.0078919081141044440000000 
-0.1037589719596667000000000 
0.0934281255397679000000000 



TABLE 4.7. NORMAL MODE AMPLITUDES FOR BIFURCATED WAVECUIDE 


p = 12 

% Noise =0.0 
Method - Optimization 


Real (x) 

-0. 6129789299973383000000000 
0. 1582655698913977000000000 
0.0297182807496905600000000 
-0. 0628349904881530000000000 
0. 0236976286523354100000000 
0. 0214179516123888800000000 
-0. 0304892422382332300000000 
0. 0060718954734892140000000 
0. 0195858898055305500000000 
-0.0200290032543383300000000 
-0.0018100034132002740000000 
0.0204861771029387600000000 
-0.0152907988556583800000000 
-0. 0089916559329415350000000 
0. 0263041570609552400000000 
-0. 0136080009193557500000000 
-0. 0260088840932548800000000 
0. 0578119165054464700000000 
-0.0164516006135241200000000 
-0. 3875515852016858000000000 
-0.7209762712782064000000000 
-0, 6777158174323345000000000 
-0. 3727708804549361000000000 
-0. 1160967435349311000000000 
-0.0160860400599178600000000 


Imaginary (x) 

-0. 7900992569759522000000000 
-0. 3230974164594634000000000 
-0. 0606695423844465700000000 
0. 1282769394906585000000000 
-0. 0483784470915562900000000 
-0.0437245114811400900000000 
0. 0622434508333472400000000 
-0. 0123957072457359000000000 
-0. 0399843769110566500000000 
0.0408889879068351700000000 
0. 0036951018453616660000000 
-0. 0418223026112559200000000 
0. 0312159955034282100000000 
0. 0183563653059957200000000 
-0. 0536996428085643900000000 
0.0277805815449464600000000 
0. 0530968450226435500000000 
-0. 1180223773132469000000000 
0. 0335857568656189800000000 
0. 7911821771790342000000000 
1.4718648726786720000000000 
1. 3835491156449770000000000 
0. 7610074828751852000000000 
0. 2370101653509874000000000 
0. 0328394647614768700000000 



TABLE 4.8. NORMAL MODE AMPLITUDES FOR BIFURCATED WAVEGUIDE 


p = 14 

% Noise =0.0 
Method - Optimization 


Real(x) 

-0. 5623470512816346000000000 
0. 1589372937710953000000000 
0.0283805678520295100000000 
-0. 0570667538152790500000000 
O. 0204014765185503600000000 
0.0173978355101469200000000 
-0. 0232337332460656100000000 
0. 0043104769692025830000000 
0. 0128448798670556800000000 
-0. 0120121195983029000000000 
-0.0009804122312131409000000 
0. 0098681612197314590000000 
-0.0064227677403004840000000 
-0. 0032103659643056290000000 
0. 0077136563489172660000000 
-0.0031235885895530720000000 
-0. 0043497089594124080000000 
0. 0062639820557906680000000 
-0. 0009198749928571885000000 
-0. 0056418959076346720000000 
0. 0061881432216453450000000 
0.0015794675372741010000000 
-0. 0190921166749397100000000 
-0. 0210777001708815600000000 
-0. 0070626586598077020000000 


Imaginary (x) 

-0. 8269013205424356000000000 
-0. 3002960644661010000000000 
-0. 0536222344740455700000000 
0.1078219036952418000000000 
-0.0385465422332256100000000 
-0. 0328714640216319800000000 
0.0438978070651003400000000 
-0.0081442136020610940000000 
-0.0242691113047228700000000 
0.0226956943610408400000000 
0.0018523905089344930000000 
-0. 0186449001859118500000000 
0.0121351750114092100000000 
0. 0060656642724994780000000 
-0. 0145741794718500400000000 
0. 0059017071345721030000000 
0. 0082183385081127160000000 
-0. 0118351653925888500000000 
0. 0017380114743673620000000 
0. 0106597960493836400000000 
-0.0116918755232731800000000 
-0.0029842453867059290000000 
0. 0360726382364162000000000 
0. 0398241989602472700000000 
0.0133441846776210500000000 



TABLE 4.9. NORMAL MODE AMPLITUDES FOR BIFURCATED WAVEGUIDE 


p = 16 

% Noise =0.0 
Method - Optimization 


Real(x) 

-0. 5446558248946477000000000 
0. 1590022040903823000000000 
0. 0279127638713525100000000 
-0. 05518X1110012006300000000 
0.0193746450703939100000000 
0.0162021929136744000000000 
-0.0211782416512421100000000 
0. 0038371934826414380000000 
0.0111371878362755800000000 
-0. 0101119924082123900000000 
-0.0007982343819241843000000 
0. 0077345572650122960000000 
-0. 0048183895891171300000000 
-0. 0022887068171741520000000 
0.0051780477628820660000000 
-0.0019508828754829160000000 
-0.0024873350411450600000000 
0.0032072018987653500000000 
-0.0004082195809085144000000 
-0. 0020634269386133310000000 
0.0017116202437072010000000 
0. 0002781881640264717000000 
-0. 0012994942829452510000000 
0. 0006190344976838535000000 
0. 0002691793983988058000000 


Imaginary (x) 

-0. 8386596642805188000000000 
-0. 2928526207490611000000000 
-0. 0514101429822789000000000 
0. 1016333894335286000000000 
-0.0356845087901091600000000 
-0. 0298414399406767200000000 
0. 0390064013395998300000000 
-0.0070674001835157070000000 - 
-0. 0205126380233020400000000 
0. 0186244183843441600000000 
0.0014701996787598670000000 
-0. 0142456220857220800000000 
0.0088745815181721930000000 
0. 0042153722634738060000000 
-0. 0095370046330763750000000 
0.0035931653836352310000000 
0. 0045812104637614720000000 
-0.0059070713608732930000000 
0. 0007518647435511409000000 
0.0038004502039422600000000 
-0. 0031524870871709220000000 
-0.0005123713449148072000000 
0. 0023934286786651200000000 
-0. 0011401471131658820000000 
-0.0004957786256640207000000 



TABLE 4.10. NORMAL MODE AMPLITUDES FOR BIFURCATED WAVEGUIDE 


p = 18 

% Noise = 0.0 
Method - Optimization 


Real(x) 

-0.5651271771500234000000000 
0.1589366664660428000000000 
0.0284631279897393900000000 
-0.0574064797458812700000000 
0. 0205920971244315600000000 
0. 0176275737885480200000000 
-0. 0236441667726119000000000 
0. 0044091006315241160000000 
0. 0132180244199123900000000 
-0. 0124499163056562900000000 
-0. 0010249649592919670000000 
0. 0104264309694853600000000 
-0. 0068763446949482210000000 
-0. 0034953938248281280000000 
0. 0085851815532825640000000 
-0. 0035811755490244880000000 
-0.0051993546230948390000000 
0.0079669691969732470000000 
-0.0012932953005069440000000 
-0.0095460782573350940000000 
0.0164964373320757300000000 
-0.0158733478044230900000000 
0. 0101404606960294600000000 
-0. 0041811699455688240000000 
0. 0008827680442359446000000 


Imaginary (x) 

-0.8250038022022767000000000 
-0. 3015211511355319000000000 
-0.0539978301244264300000000 
0. 1089067000111651000000000 
-0. 0390655785555015800000000 
-0. 0334415365484341400000000 
0. 0448557059951115800000000 
-0. 0083645714197677870000000 
-0.0250761138221265900000000 
0. 0236189243136035500000000 
0.0019444765092668660000000 
-0.0197801397123230300000000 
0.0130452174069262300000000 
0 . 0066311644325888270000000 
-0. 0162870776272356000000000 
0. 0067939022375628000000000 
0.0098637742056769140000000 
-0. 0151142576262593700000000 
0. 0024535300533226780000000 
0.0181100092808286700000000 
-0. 0312956404850722100000000 
0. 0301135679282827200000000 
-0.0192376211843818400000000 
0. 0079321606709701670000000 
-0. 0016747125932062820000000 



TABLE 4.11. NORMAL MODE AMPLITUDES FOR BIFURCATED WAVEGUIDE 

p = 10 

% Noise =5.0 
Method - Optimization 


Real (x) 

1.3490159039821090000000000 
-0. 5863026718912302000000000 
-0. 2499775949462653000000000 
0.8213025796154384000000000 
-0.4395950274506704000000000 
-0. 5356119422561171000000000 
1.0113637504109250000000000 
-0. 2590991539735201000000000 
-1.0816058668664170000000000 
1.4312553051892790000000000 
0. 1690777939542518000000000 
-2. 5417584951374730000000000 
2.6052432024487930000000000 
2.2264628193069770000000000 
-10.4485834434388700000000000 
10.7757904315912700000000000 
95. 1189043055005800000000000 
161. 6468611764827000000000000 
104.7397181855267000000000000 
-45. 4028419982236700000000000 
151.9584882904249000000000000 
144.9556870299789000000000000 
-77.0261095275353300000000000 
-23.2168301764561800000000000 
-3.1577938829878370000000000 


Imaginary (x) 

0. 1405482220186667000000000 
0.0889979696112806000000000 
0.0216084684820422300000000 
-0. 0465015532111386200000000 
0. 0199656493921559800000000 
0.0203687024467156000000000 
-0. 0343214613215909200000000 
0.0081120870832921720000000 
0.0305104364799946600000000 
-0. 0378079938575179800000000 
-0.0041949008621637160000000 
0. 0596907946134939000000000 
-0.0585530399134384800000000 
-0. 0482018569113481800000000 
0.2195731774724086000000000 
-0.2212565887805249000000000 
-1.9274887461692960000000000 
-3. 2976959960041770000000000 
-2 . 3389676240363020000000000 
0.3734639176044082000000000 
2 . 3550102529421700000000000 
2.3506298242204450000000000 
1.2635392214914600000000000 
0. 38244948326031590.00000000 
0.0521985610648173600000000 



TABLE 4.12. NORMAL MODE AMPLITUDES FOR BIFURCATED WAVEGUIDE 


p - 12 

% Noise =5.0 
Method - Optimization 


Real(x) 

-1.0023585290789630000000000 
0.3100646966322169000000000 
0.0764043293211506000000000 
-0.2188818087689583000000000 
0. 1032826582243121000000000 
0. 1104462254567384000000000 
-0. 1893837195006025000000000 
0.0430740826754165400000000 
0.1591954780580723000000000 
-0.1891541368826273000000000 
-0.0189840947476063100000000 
0.2490634533462135000000000 
-0.2101299921672098000000000 

-0.1407034729495921000000000 
0.4612647131903209000000000 
-0.2706531744915136000000000 
-0.5803486017311115000000000 
1. 4598242886131840000000000 
-0.4709115339550887000000000 
12.6652815359736800000000000 
27.0186318353265900000000000 
29 . 3480633346051800000000000 
18.8453451106995300000000000 
-6 .9492696705619030000000000 
-1.1624752439276550000000000 


Imaginary (x) 

-0.9784694148900999000000000 
-0 . 3903625019492461000000000 
-0.0732690924796732000000000 
0 . 1566603139235494000000000 
-0.0593859156631996500000000 
-0.0527159655088688100000000 
0 .0749584454092393000000000 
-0.0149428106380085500000000 
-0 .0475537426564181200000000 
0.0482202879597424600000000 
0.0042728301146030010000000 
-0.0485423423836005200000000 
0.0359214787383149000000000 

0 . 0209361149793653300000000 
-0.0606930921042045200000000 
0.0310092901230144100000000 
0.0584520098722679400000000 
-0 . 128192566 729 7322000000000 
0.0359592539754477200000000 
0 . 8307820610111580000000000 
1 .5152024180475260000000000 
1.3936116065156360000000000 
0 . 7489315044303903000000000 
0.2280872481229062000000000 
0.0311308187088503600000000 
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TABLE 4.13. NORMAL MODE AMPLITUDES FOR BIFURCATED WAVEGUIDE 
p = 14 

% Noise = 5.0 
Method - Optimization 


Real (x) 


Imaginary (x) 


-0.5535550281324885000000000 

0.1635719855462542000000000 

-0 . 829 544634 742 94 94000000000 
-0.3038167382678901000000000 

0.0311988610014427800000000 

-0.0611417713221645700000000 

-0.0532273252857062000000000 

0.1076639475740755000000000 

0.0 214448255835 784900000000 
0.0228339394367208700000000 

-0.0392169717302455000000000 

-0.0329356018593202500000000 

-0.028931297003742 8600000000 
0.0046510908638282310000000 

0.0437287867613229900000000 

-0.0080222310723188200000000 

0.015 293433 681 5 745400000000 
-0.0134289209924681100000000 

-0.0239326092889081100000000 : 

0.022 6464692343 5 1 01 00000000 

-0.0011709794127824960000000 

0.0115087258743531200000000 

0.0019512956885194050000000 

-0.0184954540245556200000000 

-0.0075709684507100940000000 

0.0 1 2 1 3 3 871 8 5 29 5 6 74 0 0 000 0 00 

f 

-0*003320647488421 7900000000 
0.0090045487049005090000000 

0.0059410557979731360000000 

-0.0144562914763565400000000 

-0.0025603285071671670000000 
-0.005521742497811 8710000000 

0.0058493771931077940000000 

0.0081905055584625310000000 

0.0093050164375639800000000 
-0.001633583 6529978510000000 

-0.0116898921998918000000000 j 
0.0017496585025476980000000 

-0.0085204231626269930000000 

0.0100238654750026300000000 

0.0105138293487714900000000 1 

-0.0114963364050298100000000 

0.0036803951979292240000000 

-0.0390000126188903100000000 

-0.002921507777592 5260000000 
0.0352751703289899600000000 

-0.0493315492269055200000000 
-0. 021 3 33 8 Q4 509955 75 00000000 

0.0387548632824920100000000 

0.0128830485410871100000000 
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TABLE 4.14. NORMAL MODE AMPLITUDES FOR BIFURCATED WAVEGUIDE 

p = 16 

% Noise = 5.0 
Method - Optimization 


Real (x) 

Imaginary (x) 

-0.5331505330900437000000000 
0.1593085632503 803000000000 

-0.8524684163381 835000000000 
-0.2991039010335467000000000 

0.0340763500362040100000000 

-0.0511506396512322600000000 

-0.05422018745991472 00000000 
0.1043751398286410000000000 

0.0188311876322114200000000 

0.0178151606320703800000000 

-0.0364351932761602000000000 
-0.030 201 85652 02011100 0000 00 

-0.0228754907403309600000000 

0.0012913377923405920000000 

0.0394376387480535400000000 

-0.0074105211093037410000000 

0.0116050329218668500000000 

-0.0107950992904147800000000 

-0.02077922447552 90300000000 
0.0188600312642632500000000 

-0.0004262007793781324000000 

0.0101703207249470800000000 

0.0013685322713389020000000 

-0.0143703042285351400000000 

-0.004764187654901 7790000000 

0.0089254542739283100000000 

-0.00 13812417*71 4357780000000 
0.0067469808732455280000000 

0.0042829158099681870000000 
-0.00962 55926622258750000000 

-0.00 1 65 1 374 6103 943110000000 
-0.0022281887242793660000000 

0.003619194245966 7690000000 
0.0046635431166843470000000 

0.0034818223541 944810000000 
-0.0009351837287743341000000 

-0.0059715622805034380000000 

0.0007733212388813635000000 

-0.0015920753768636370000000 

0.0011454630165291140000000 

0.0038168633330529200000000 

-0.0031828482580409380000000 

0.0012829285790696270000000 

-0.0013879393604450130000000 

-0.0005586440510891364000000 

0.0024187133443262950000000 


0*0001280284621 738308 000000 -0.0011558709048975160000000 

-0.0002602833834927673000000 -0.0004996919010203439000000 
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TABLE 4.15. NORMAL MODE AMPLITUDES FOR BIFURCATED WAVEGUIDE 

p = 18 

% Noise = 5.0 
Method - Optimization 


Real(x) 


Imaginary (x) 


-0.5756673958387165000000000 
0.152 8958902455 5180 00000000 

-0 . 844600865*4762 2 31000000000 
-0.3121877427107791000000000 

0.0281476301297940200000000 

-0.0538862103246514800000000 

-0.0574780132759763000000000 

0.1113921309626813000000000 

0.0203037128513620800000000 

0.0176468163534723300000000 

-0.0393079969447663000000000 

-0.0342918825747588100000000 

-0.0264026722677921400000000 

0.0064070701861388890000000 

0.0455691590856567000000000 

-0.0083140685090928410000000 

0.0147760271128326600000000 

-0.0122956545024289200000000 

-0.0256588566117524200000000 

0.0241537726341783200000000 

-0.0005083578254922163000000 

0.0114715450530849300000000 

0. 001 9421 8B1 7191 465 80000000 
-0.0201286945337159600000000 

-0.0091135931936450570000000 

0.0132990588935517800000000 

- 0 . 0 041 3776 86 72 03 06 9 80 00000 0 
0.0084795681701519340000000 

0.0066388490766845080000000 

-0.0165578391221838600000000 

-0.0045609385308836570000000 

-0.0062437489536120490000000 

0.0068705984302362800000000 

0.0100731914298415500000000 

0.0104395129488849200000000 

-0.0025787959228411210000000 

-0.0154539361167026100000000 

0.0025173604248903590000000 

-0.01240193 04870910400000000 
0.02220931275452491 00000000 

0. 0185906849200285900000000 
-0.0321064875567949100000000 

-0.0222559640900727100000000 

0.0153254116355631300000000 

0.0309185110096987100000000 
-0.01979087308641 07000000000 

-0.0078002730530776860000000 

0.0026957561223150330000000 

0.0081467347535768300000000 
-0.00172 302911 72 084940000000 
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5 . SECOND EXAMPLE 

In the next example, an attempt was made to solve an integral 

4 

equation related to an electromagnetic remote sensing problem. The 
equation is 


/ f(x) sinh 2 y(x - 1) dx = g(y) . (5) 

o 


The method of moments was used to convert this continuous operator 
equation to a matrix equation. 

Let f(x) be approximated by a sum of expansion functions, so that 


N 

f(x) - l a <j» (x) (6) 

j-1 3 3 


where <)>j (x) is an expansion function and a_. is its respective coefficient. 
We note that the operator of Equation (5) is linear which permits us to 
distribute the integration over the series approximating f(x). Thus, 


’ 2 

l a. / <h(x) sinh y(x - 1) dx » g(y') . 

j=l 3 o 3 

We next take the inner product of both sides with the testing functions 
6(y - y^) which results in 

N 1 

l a, I <K(x) .sinh y, (x - 1) dx = g(y.) . (7) 

j=l 3 o 3 


To finally get this into matrix form, we define 
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A ij = ^ sinh 2 y i (x - 1) dx 


gCy i ) = 6 i 


Equation (7) can now be written as 

E 


T a.A. . = B 
3-1 3 13 1 


Expressed in matrix form, we get 


- ; 





A 


a 

a 

B 



i 




To effect a solution on the computer, we must specify the expansion 

functions . (x) . In this paper, 

1 

fjOO = 


A_, . is then found to be 
ij 


Aj. = / x j 1 : sinh 2 y (x - 1) dx 

J A 


To carry out the integration, we make the substitutions 


u = y ± (x - 1) 
du = y. dx 

l 


A. . 



(— + d ^" 1 


y i 


sinh^ 


u du 
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Using the binomial expansion theorem, 


A . = / 


~y± ' x 


(i - l)t fu_ ,k 


k t Q (j - k - l)!k! \y ± 


sinh^ u du 


° ^ 
J y . 

- y i * 


i , . , v . k-l 

y r-JLIi 

(j - k)t (k - 1)! \y ± 


sinh u du 


. f v ~ k 0 - HI S(k) 

k £ 1 y i (j - k) ! (k - 1> ! S(k) 


where 


S(k) = / x k 1 sinh Z x dx 


-y* 


From an integral table we find 


/ x n sinh 2 (x) dx 


-x 


n+1 




2 (n + 1) 


+ 4 I 


k=o 


x n 2k ainh (2x) 


2 2k (n - 2k)! 


n-2k-l f . 

: oosh(2x) 

2k+l * 

2 ZK ^(n - 2k - 1)! 


We now have an explicit formula for the A matrix* 

Before going further, we must choose the domain of g(y), the number 
of expansion functions, and the number of testing functions. Although 
it is possible to have more testing functions than expansion functions, 
it is simpler to generate the solution if there are an equal number of 
the two types of functions and, thus, a square matrix. In this 
section, solutions were sought for varying numbers of expansion functions 
and for varying domains of g(y). 
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The testing procedure was as follows. An f(x) was assumed from 
which the right-hand side of Equation (5) could be calculated exactly. 
The inner product of this g(y) was then taken with the testing functions 
< 5 (y _ y^) to find the elements of the matrix. Then it was attempted to 
reconstruct the assumed f(x) given the A and 3 matrices. A solution 

2 

was attempted for each of four assumed functions: f(x) = 1, f(x) = x , 

f(x) = sin (ttx), f(x) = sin (5irx) . 

The coefficients were found by using both Gaussian elimination and 
optimization techniques. The Gaussian elimination method resulted in 
the more exact solution. The same problems of overflow were again 
encountered when using DFMFP for larger domains of g(y)« When the 
function value and gradient were multiplied by a number within FUNCT, 
answers were obtainable for larger domains of g(y)» but they were much 
less accurate. The answers reported were found without using any 
multiplying factor. 

After the coefficients for the expansion functions were found, the 
power series approximation to f(x) was evaluated at various values of x. 
These values along with the true values are tabulated in the following 


tables. 



I 

I 
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TABLE 5.1. RECONSTRUCTED FUNCTION VALUES 

Domain of g(y) is { y : 0.0 < y < YE} 

N * 4 

Method - Optimization 


YE= 2.50 


CORRECT 

1.00000000000 

1.00000000000 

1.00000000000 

1.00000000000 


CORRECT 

0.382683359860 

0.923879442413 

0.923879682675 

0.382683939902 


CORRECT 

0.156250000000D-01 

0.140625000000 

0.390625000000 

0.765625000000 


F (X)=l 

CALCULATED 

1.00032439684 

0.999317671743 

1.00011673488 

1.00894447917 

F(X)=SIN( 

CALCULATED 

0.343217595064 

0.749150436478 

1.33387338504 

2.15369604060 

F(X)=X**2 

CALCULATED 

0.205291063448D-01 

0.143298227823 

0.372290024440 

0.747747362151 


ERROR 

-0. 32439684 0952D-03 
0.682328256995D-03 
-0.116734879538D-03 
-0.894447916742D-02 

PI*X) 

ERROR 

■ 0.394657647956D-01 
0.174729005936 
-0.409993702365 
-1.77101210070 


ERROR 

-0.4 904106344 7 7D-02 
-0. 267322782335D-02 
0.183349755601D-01 
0.178776378488D-01 


F(X)=SIN (5*PI*X) 


CORRECT 

0.923879682675 

-0.382684519944 

-0.382681619733 

0.923878481365 


CALCULATED 
0.395055487014 
0. 181708889698D-01 
0.486904539040 
-1.01009604106 


ERROR ' 

0.528824195661 

-0.364513630974 

0.104222919307 

1.93397452243 



TABLE 5. 


YE= 5 . 00 


CORRECT 

1.00000000000 

1.00000000000 

1.00000000000 

1.00000000000 


CORRECT 

0.382683359860 

0.923879442413 

0.923879682675 

0.382683939902 


CORRECT 

0. 156250000000D-01 
0.140625000000 
0.390625000000 
0.765625000000 


CORRECT 

0.923879682675 

-0.382684519944 

-0.382681619733 

0.923878481365 


2. RECONSTRUCTED FUNCTION VALUES 

Domain of g(y) is {y: 0.0 < y < YE} 

N = 4 

Method - Optimization 


F(X)=1 

CALCULATED 

1.00036753360 

0.999159300662 

1.00044495943 

1.01203072637 

F (X)=SIN( : 

CALCULATED 
0.386640347517 
0.916829181261 
0. 946378140612 
0.286561483641 

F(X)=X**2 

CALCULATED 

0. 181273364432D-01 

0.134899954776 

0.393660842679 

0.847570483165 


ERROR 

-0. 367533598703D-03 
0.840699338224D-03 
-0. 444959430271D-03 
-0. 120307263687D-01 

*X) 

ERROR 

-0. 395698765678D-02. 

0. 705026115271D-02 
-0. 224984579378D-01 
0.961224562607D-01 


ERROR 

-0.250233644318D-02 
0.572504522359D-02 
-0. 303584267861D-02 
-0. 819454831652D-01 


F(X)=SIN(5*PI*X) 


CALCULATED 

0.383414934051 

-0.112563400807 

-0.380118460258 

-0.314911510689 


ERROR 

0.540464748623 

-0.270121119137 

-0.256315947507D-02 

1.23878999205 



TABLE 5. 


YE= 10.0 


CORRECT 

1.00000000000 

1.00000000000 

1.00000000000 

1.00000000000 


CORRECT 

0.382683359860 

0.923879442413 

0.923879682675 

0.382683939902 


CORRECT 

0. 156250000000D-01 
0.140625000000 
0.390625000000 
0.765625000000 


CORRECT 

0.923879682675 

-0.382684519944 

-0.382681619733 

0.923878481365 


3. RECONSTRUCTED FUNCTION VALUES 

Domain of g(y) is {y: 0.0 < y < YE} 

N - 4 

Method - Optimization 


F(X)=1 

CALCULATED 

1.00014282939 

0.999126782873 

1.00849210828 

1.04374297373 

F (X)=SIN ( 

CALCULATED 

0.384247884274 

0.917885451900 

1.02530110059 

0.586273513597 

F(X)=X**2 

CALCULATED 

0.161546194444D-01 

0.137374034664 

0.421627555424 

0.925855142525 


ERROR 

-0. 14 28293893 66D-03 
0.87 3217127 032D-03 
-0.84 9210827 702D-02 
-0.437429737327D-01 

PI*X) 

ERROR 

-0.156452441388D-02 
0. 599399051295D-02 
-0.101421417919 
-0.203589573695 


ERROR 

-0. 529619444390D-03 
0, 325096533635D-02 
-0. 310025554243D-01 
-0.160230142525 


F(X)=SIN(5*P I*X) 


CALCULATED 

0.699458232406 

0.593549636234 

-5.43017826864 

-14.7377692041 


ERROR 

0.224421450269 

0.210865116290 

5.04749664890 

15.6616476855 
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TABLE 5.4. RECONSTRUCTED FUNCTION VALUES 

Domain of g(y) is {y: 0.0 < y < YE} 

N = 8 

Method - Optimization 

YE* . 2.50 __ 


F1XJ=_L_._ 

CORR ECT 

1. 00000000000 
_1. OOQOQOOO_OOQ 

1.00000000000 
1.00000000000 
1.00000000000 
1 .00000000000 
1.00000000000 
1.00000000000 


F( X)«S IN< PI»X) 

CALCULAJEO ERROR 

0.27342 2296279 -0.7833201274840-01 

0.420492266405 _ 0.13507 7868735 

0.602170243084 0.229299260217 

0.837269729007 __ 0.14351549 7 809 

1.15972775755 -0.178942408249 

1 .,62 97 8 56056 0 -0 .79 331 5753496 

2.35054786118 -1.79497720402 

3, 4910794474 7 . -3.2953335A&11 


F <X)*X*»2 


CORRECT 


CALCULATED 

ERROR 

0.3 9062 5000000D 
0.3515625000000 

-02 

-oj^ 

0.1 14546027231D-01 
0.3379683235590-01 

-0.7 548352723090-02 
0. 1359417644150-02 

0.9765625000000 

0.191406250000 

-01 

0.7554476107380-01 
0. 152341162766 

0.2211 14889262 D-Ol 
0.3 9065 087234 ID- 01 

0.316406250000 

0.472656250000 


0.293917244040 

0.554855542663 

0.2248900595950-01 
-0.821 992926629D-01 

0.660156250000 

0.878906250000 


1.03058816351 

1.87976442315 

-0.370431913514 

-1.00085317315 



FIX) 

*SIN( 5*PI*X1 

CORRECT 


CALCULATED 

ERROR 

0.831469503301 

0.195090899300 


0.505532915032 

0.275134140227 

0.325936587469 

-0.8004324092670-01 

-0.980785471784 

0.555569091092 


0.4257988 46 048D-01 
-0.1777649^9109 

-1.02336535639 

0.733334040200 

0.555571701211 

-0.980784859362 


-6. 356365250965 
-0.436689979391 

0.911936952176 

-0.544094879970 

0.19508782 0453 
0.831 471247324 


-0.317627525943 

0.171103683967 

0.512715346396 

0.660367563357 


OR REC T 

0.195090283531 

0.555570135140 

0.831469503301 

0.980785226816 

0.980785349300 

0.831 469852 106_ 

0.555570657164 

0.195090899300 


CAL CULATED 


ERROR 


1.00393485659 

Q ,99 7 453.8 11 56 2 

0.988221451970 

0.9815 306576 63 

0.990776725528 

1.04472094405 


-0.3934856589550-02 


0.1177854803040-01 
0. 1 84 6934233740-01 
0. 92232 74471 90D-02 
-0.4472 Q9440476D-01 


1.19851648324 

1.54933565809 


-0.198516483237 
-O t 549 335658094 



36 


TABLE 5.5. RECONSTRUCTED FUNCTION VALUES 

Domain of g(y) is {y: 0.0 < y < YE} 

N = 8 


5.00 


Method - Optimization 


F ( X )=1 


CORRECT 

1 . 00000000000 

i. oooooo ooooa 

1. 00000000000 

1.00000000000 

1 . 00000000000 

1.00000000000 

1 . 00000000000 

1.00000000000 


CALCULATED 
1.00072335365 
. .1,002 85 90.5 .3 *5 0 
0.996956056567 
0,967057559989 
0.985204565105 
1.01880960839 
1. 16191829994 
1.45117868134 


ERROR ; 

-0.723353645402D-03 
- 0 . 2.85 9 05 3 899230*02 
0. 3 04 5 943 4 3 34 8 D- 02 
... Q, 1 29424400106D -0 1 
0.1 4795454 8953 D-01 
-0.18809 60 839020-01 
-0.141918299943 
- 0, 4 51178681338 


F(X1=S1N( PI*XI 


CORRECT 

0.195090233531 
0.555570135140 
0.931469503301 
0.980785226816 
0.980785349300 
0.831469852106 
0.555570657164 
0.1 95090899300 


CALCULATED 
0.197705476216 
0.554010517714 
0.825475048767 
0.985602069239 
1.00262 877598 
0.862362751656 
0.47! 787692286 
■0.1 35632 573Z67 


ERROR 

-0.26151 92685510-02 
_0 , 1_55 9 617 42 567D-02 
0. 5994454553640-02 
-0.6616822 62 309D-02 
-0 . 2 166362667760-01 
- 0 .1 08 72 89 93 498 D-01 
0. 8378296487980-01 
0.330723472567 


C ORRECT 

0. 39062 50000000- 02 
-Q- «351 5625Q0000D-01 
0.976562500000D-01 
-Qj. 1 914062 50000 
0.316406250000 
0.472656250000 
0.660156250000 
0.878906250000 


CALCULATED 

0. 58791 4934171D-02 
0 .4071 65070652 D-01 
0. 89872784442 7D-01 
0.165785868885 
0.292966269075 

0.51754757 2 02 0 

0.921556305699 

1.64293166159 


)=X*»2 

ERROR 

! -0. 1 9728 99 34171 D- 02 
L -0.5560257 06516 D-Q2 
L 0. 7 7B 3 46 5 5 573 3D- 02 
0.256203 81 1.1 5 OD^O 1 
0.2343998092470-01 
-0.44891 3 2202 0 00- Q1 
-0.261400055699 
-0.764025411595 


FCX)»SIN(5*PI*X> 


CORRECT 

0.831469503301 
0.1 95090899300 
■0.980785471784 
0.555569091 092 
0.555571701211 
■0.980784859367 
0.195087820453 
0.831471247324 


C ALCULATED 

9.540639791.880 
0.309035825687 
■0.292765693744D-01 
*0. 3 96834904775 
•0.622275694373 
^Q_. 3 6998 1297206 
0.963790692757 
4. 39311430002 


ERR OR 

0.290829711421 
3944926387 
*0.951 50890241 0 
0.95240399 5 866 
l. 17784739558 
*0.61090356 2156 
■0.768702872305 
-3.56164305270 
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TABLE 5.6. RECONSTRUCTED FUNCTION VALUES 

Domain of g(y) is {y: 0.0 < y < YE} 

N * 8 


10.0 


Method - Optimization 




CORRECT - 

1. 00000000000 

1 . 00000000000 

1.00000000000 
1.00000000000 
1.00000000000 
1.000000 000 00 

1 . 00000000000 

1.00000000000 


C A.LCULATE.D 

1.00005886241 
0.999973701066 
0.998955318116 
0.998873252508 . 

1.00798604263 
1.0 475 5181 393 
1. 16201968226 
1. 43361 575539 ... 


ERROR 

-0.5 8362 4053408D-04 
0. 2 62 9 8 9335540D-04 
0. 1044681883770-02 
0.112 6 74749162D-02 
-0.7986042 6341 5D- 02 
-0.4755181393130-01 
-0.162019682260 
-0.4336157 55386 


F ( X 1= S I N ( PI*X) 


c 0 R R EC T CALCUL AT ED ERR3R 


0.195090283531 

0.555570135140 

0.831469503301 

0.980785226816 

0.980785349300 

0.831469852106 

0.1 9 6l 3 89401 03 
0.553063976003 
0.833753063662 
0.998249004789 
0.991528839861 
0.73891632551 6 

-0.1048656572440-02 
0. 2506159136900-02 
-0. 2283560360790-02 
-0.1746377797230-01 
-0.1074349056030-01 
0. 9255352658970-01 

0.555570657164 

Q.195090899300 

0.140665754719 

-0.9343.18762295 

0.414904902445 
1.U 2-9-4D966159 


f <X) = 

X**2 

CORRECT 

CALCULATED 

ERROR 

i 0.3906250000000-02 
0.351 5 62 5 OOOOOD-Ol 

0. 3869318056310-02 
0. 352509547590D-01 

0.3693194368870-04 
-0 • 9 47 04 75 90475 D- 04 

0.9765625000000-01 
: 0.191406250000 

0.9 7687 8746 8 OOD- 01 
0.190475302882 

- 0. 3 162 46 8 0022 4 D- 04 
0.9309471178710-03 

0.316406250000 

0.472656250000 

0.314587700686 

0.476087901865 

0.1 81854931389D-02 
-0. 343165186505 D-02 

0.660156250000 

0.878906250000 

0.691732792335 

0.997311709271 

-0.3157654233540-01 
-0.1 18405459271 


FIX) 

»SIN(5*PI*X) 

CORRECT 

CALCULATED 

ERROR 

0.831469503301 

0.195090899300 

0.843253556418 

0.276768128120 

-0.1178405311690-01 
-0.8 16 772 2 882 OOD- 01 

-0.980785471784 
0.555569091 Q92 

-0.996401167874 

-0.417389864164 

0.1 56156960898D-01 
0.972958955255 

0.555571701211 
-0.980784859362 

2.55103363604 
3.68202472L97 

-1.99 546i93483 
-4.66280958134 


-13 .2533522909 
-67.5532731 7^6 

10.4484401114 

68.3847444209 


TABLE 


CORRECT 

1.00000000000 

1.00000000000 

1.00000000000 

1.00000000000 


CORRECT 

0.382683359860 

0.923879442413 

0.923879682675 

0.382683939902 


CORRECT 

0. 156250000000D- 
0.140625000000 
0.390625000000 
0.765625000000 


CORRECT 

0.923879682675 

-0.382684519944 

-0.382681619733 

0.923878481365 


.7. RECONSTRUCTED FUNCTION VALUES 

Domain of g(y) is { y : 0.0 < y < YE} 

N = 4 


Method - Gauss 
YE = 2.5 


F(X)=1 


CALCULATED 

1.00000000000 

1.00000000000 

1.00000000000 

0.999999999992 


ERROR 

0. 444 699832 5 14D-12 
-0. 591304782915D-12 
-0. 352334788016D-12 
0. 830882584957D-11 


F(X)=SIN( P*X) 


CALCULATED 

0.386530753386 

0.918154180182 

0.943836456584 

0.266374767586 

F(X)=X**2 

CALCULATED 

01 0. 15 624 999999 6D-01 

0.140625000001 
0.390625000000 
0.765624999992 


ERROR 

-0.384739352603D-02 
0.572526223089D-02 
-0. 199567739099D-01 
0.116309172316 


ERROR 

0. 429978101735D-12 
-0. 570279934387D-12 
-0. 344529960117D-12 
0.800945421098D-11 


F(X)=SIN (5*PI*X) 


CALCULATED 

0.412510282984 

-0.273710561743 

-0.132218047460 

2.10208631203 


ERROR 

0.511369399691 

-0.108973958200 

-0.250463572273 

-1.17820783067 
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TABLE 5.8. RECONSTRUCTION FUNCTION VALUES 


Domain of g(y) is {y: 0.0 < y < YE} 

N = 4 


Method - Gauss 
YE = 5.0 


CORRECT 

1.00000000000 

1.00000000000 

1.00000000000 

1.00000000000 


CORRECT 

0.382683359860 

0.923879442413 

0.923879682675 

0.382683939902 


CORRECT 

0. 156250000000D-01 
0.140625000000 
0.390625000000 
0.765625000000 


CORRECT 

0.923879682675 

-0.382684519944 

-0.382681619733 

0.923878481365 


F(X)=1 

CALCULATED 

1.00000000000 

1.00000000000 

1.00000000000 

1.00000000000 

F(X)=SIN( PI*X) 

CALCULATED 

0.385155669915 

0.920226522169 

0.944574260642 

0.237926744807 

F(X)=**2 

CALCULATED 
0 . 156250000000D-01 
0.140625000000 
0.390625000000 
0.765625000000 

F (X) =SIN (5*PI*X) 

CALCULATED 

0.5346S4591614 

-0.458650735253 

-0.196595332038 

4.63881812009 


ERROR 

0.525968157916D-14 
-0 . 910382880193D-14 
-0 . 155431223448D-14 
0 . 119002030452D-12 


ERROR 

-0 . 247321005537D-02 
0.365292024469D-02 
-0. 206945 7 796 74D-01 
0.144757195095 


ERROR 

0 . 95 149 5 8265 7 3D- 15 
0.159594559790D-14 
-0.661970478433D-14 
-0, 2 70339306 496D- 13 


ERROR 

0.389195091061 
0. 759662153098D-01 
-0.186086287695 
-3.71493963872 



TABLE 5.9. RECONSTRUCTED FUNCTION VALUES 


Domain of g(y) is {y: 0.0 < y < YE} 

N = 4 


Method - Gauss 
YE = 10.0 


CORRECT 

1.00000000000 

1.00000000000 

1.00000000000 

1.00000000000 


CORRECT 

0.382683359860 

0.923879442413 

0.923879682675 

0.382683939902 


CORRECT 

0.156250000000D-01 

0.140625000000 

0.390625000000 

0.765625000000 


CORRECT 

0.923879682675 

-0.382684519944 

-0.382681619733 

0.923878481365 


F(X) = 1 

CALCULATED 

1.00000000000 

1.00000000000 

1.00000000000 

1.00000000000 

F(X)=SIN ( PI*X) 

CALCULATED 

0.382767600379 

0.926951138269 

0.937890041604 

0.135359441899 

F(X) =X**2 

CALCULATED 

0.156250000000D-01 

0.140625000000 

0.390625000000 

0.765625000000 

F(X)=SIN (5*PI*X) 

CALCULATED 

0.801707075221 

-1.22038659825 

0.579689758030 

16.2945681660 


ERROR 

0 . 5 8286 708792 8D-15 
-0. 177635683940D-14 
0.459354776439D-14 
0 . 367206265395D-13 


ERROR 

-0. 842405191514D-04 
-0.307169585562D-02 
-0.140103589297D-01 
0.247324498003 


ERROR 

-0.160548657702D-14 
0. 846545056277D-15 
O.945077349712D-14 
0 . 14 16922 135 18D-13 


ERROR 

0.122172607453 
0. 837702078308 
-0.962371377763 
-15.370 6 89 6 846 
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TABLE 5 


CORRECT 

1.00000000000 

1.00000000000 

1.00000000000 

1.00000000000 

1.00000000000 

1.00000000000 

1.00000000000 

1.00000000000 


CORRECT 

0.195090283531 

0.555570135140 

0.831469503301 

0.980785226816 

0.980785349300 

0.831469852106 

0.555570657164 

0.195090899300 


CORRECT 

0. 390625000000D-02 
0. 351562500000D-01 
0. 976562500000 D-Ol 
0.191406250000 
0.316406250000 
0.472656250000 
0.660156250000 
0.878906250000 


CORRECT 

0.831469503301 
0.395090899300 
-0.980785471784 
0,555569091092 
0.555571701211 
-0.980784859362 
0. 195087820453 
0.831471247324 


.10. RECONSTRUCTED FUNCTION 
Domain of g(y) is {y*. 

N = 8 

Method - Gauss 
YE = 2.5 


F(X) = 1 

CALCULATED 

1.00000013027 

0.999999857246 

1.00000024124 

0.999999892182 

0.999999411834 

1.00000110476 

1.00000158535 

0.999978527908 

F(X)=SIN( P1*X) 

CALCULATED 

0.195085993822 

0.555577915809 

0.831459612337 

0.980779141549 

0.980836579725 

0.831411302485 

0.555367676864 

0.197091066722 

F(X)=**2 

CALCULATED 

O.390586223766D-02 

0.351566739390D-01 

0.976555334236D-01 

0.191406572877 

0.316407986865 

0.472652977554 

0.660151602768 

0.878969452133 

F(X)=SIN(5*PI*X) 

CALCULATED 

0.898234618366 

0.149974476070 

-0.839473802652 

0.173525822960 

1.14585503368 

-1.06724381390 

-2.54522335865 

19.8919035871 


VALUES 

0.0 < y < YE} 


ERROR 
-0.130265724207D-06 
0.142754472132D-06 
-0.241237670728D-06 
0.107818383319D-06 
0.588165838311D-06 
-0 . 11047636 76 32D-0 5 
-0.158534871408D-05 
0.214720919098D-04 


ERROR 

0.428970823621D-05 
-0.778066926425D-05 
0.989096349792D-05 
0 .60852676 79 40D-05 
-0 . 512 3042 46 149D-0 4 
0 . 5 85496207 7 40D-04 
0.202980300447D-03 
-0.200016742257D-02 


ERROR 

0.387762337962D-06 
-0.423938958122D-06 
0.716576379609D-06 
-O.322876702860D-06 
-0 . 17 36 8650016 7D-05 
0. 327244558852D-05 
0.464723234163D-05 
-0.632021333217D-04 


ERROR 

-0.667651150647D-01 
0.4511642 32299D-01 
-0.141311669133 
0.382043268131 
-0.590283332468 
0.864589545368D-01 
2.74031117911 
-19.0604323397 


TABLE 5.11. RECONSTRUCTED FUNCTION 

Domain of g(y) is {y: 0.0 < y < YE} 

N = 8 

Method - Gauss 


CORRECT 

1.00000000000 

1.00000000000 

1.00000000000 

1.00000000000 

1.00000000000 

1.00000000000 

1.00000000000 

1.00000000000 


CORRECT 

0.195090283531 

0.555570135140 

0.831469503301 

0.980785226816 

0.980785349300 

0.831469852106 

0.555570657164 

0.195090899300 


CORRECT 

0. 390625000000D-02 

0.351562500000D-01 

0.9 76562 500000D-01 

0.191406250000 

0.316406250000 

0.472656250000 

0.660156250000 

0.878906250000 


CORRECT 

0.831469503301 
0. 195090899300 
-0.980785471784 
0.555569091092 
0.555571701211 
-0.980784859362 
0.195087820453 
0.831471247324 


YE = 5.0 

F (X)=l 

CALCULATED 

0.999999999284 

0.999999999250 

0.999999999488 

1.00000000650 

0.999999984727 

1.00000000406 

1.00000009559 

0.999999241948 

F(X)=SIN ( PI*X) 

CALCULATED 

0.195089881285 

0.555570264242 

0.831468443880 

0.980789305854 

0.980778477412 

0.831464966588 

0.555644633901 

0.194541305072 

F (X)=X**2 

CALCULATED 
0 . 390625207842D-02 
0. 351562513982D-01 
0.9 7656252 133 9D-01 
0.191406233669 
0.316406284462 
0.472656245736 
0.660156026616 
0.878907968205 

F (X) =SIN(5*PIX) 

CALCULATED 
0.832577376712 
0. 148003583170 
-0.927927000837 
0.524168699037 
0.562384312986 
-1.19389945458 
1.65577018395 
-10.6135982335 


ERROR 

0.715634690396D-09 
0. 750 4 85 85 40 1 7D-09 
0.5122 441082 42D-09 
-0.65016689 8097D-08 
0.152729624225D-07 
-0.406042865997D-08 
-0.95586 362 380 8D-0 7 
0.758052063662D-06 


ERROR 

0. 402245829 756D-06 
-0. 129101766139D-06 
0. 105942092377D-05 
-0.407903726130D-05 
0.687188812204D-05 
0. 488551805439D-05 
-0. 739767369056D-04 
0.549594228093D-03 


ERROR 

-0.207841641554D-08 
-0. 139816589903D-08 
-0.213385389913D-08 
0.16 33100 80 87 5D-0 7 
-0. 344617861608D-07 
0.426397697706D-08 
0.223383635645D-06 
-0.171820503503D-05 


ERROR 

-0.110787341147D-02 
0. 470873161298D-01 
-0 .528584709 46 8D-01 
0.314003920541D-01 
-0. 6812611774930-02 
0.213114595217 
-1.46068236350 
11. 4450694808 
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TABLE 5.12. RECONSTRUCTED FUNCTION VALUES 

Domain of g(y) is {y: 0.0 < y < YE} 


CORRECT 

1.00000000000 

1.00000000000 

1.00000000000 

1.00000000000 

1.00000000000 

1.00000000000 

1.00000000000 

1.00000000000 


CORRECT 

0.195090283531 

0.555570135140 

0.831469503301 

0.980785226816 

0.980785349300 

0.831469852106 

0.555570657164 

0.195090899300 


CORRECT 

0.390625000000D-02 

0. 351562500000D-01 

0.976562500000D-01 

0.191406250000 

0.316406250000 

0.472656250000 

0.660156250000 

0.878906250000 


CORRECT 

0.831469503301 

0.195090899300 

-0.980785471784 

0.555569091092 

0.555571701211 

-0.980784859362 

0.195087820453 

0.831471247324 


N = 8 

Method - Gauss 

YE ■ 10.0 

F(X)=1 

CALCULATED 

1.00000000000 

1.00000000000 

0.999999999994 

0.999999999988 

1.00000000007 

0.999999999925 

0.999999999631 

1.00000000363 

F(X)=SIN( P*X) 

CALCULATED 

0.195090186347 

0.555570455895 

0.831468790285 

0.980786902937 

0.980783300617 

0.831464802191 

0.555613007649 

0.194778883599 

F(X) =X**2 

CALCULATED 
0 . 390625000018D-02 
0 . 351562499950D-01 
0.9 765625 000 71D-01 
0.191406250009 
0.316406249930 
0.472656250094 
0.660156250365 
0.878906246126 

F(X)=SIN(5*PI*X) 

CALCULATED 

0.844347164257 

0.164836718532 

-0.924952706385 

0.406505386032 

0.875356958384 

-1.32955040482 

-0.244829169278 

5.11032578941 


ERROR 

-0.945910016981D-13 
-0. 49624 74 87547D-11 
0.585063941738D-11 
0.124843885230D-10 
-0 . 6 90 3 89 9 6 7 4 17D-10 
0. 750776524283D-10 
0.3689 19728100D-09 
-0. 362715968372D-08 


ERROR 

0.971837593733D-07 
-0.320754682342D-06 
0.713015537293D-06 
-0.16 7612 106 55 5D-05 
0.204868387196D-05 
0.50499 1541872D-05 
-0. 42350485029 7D-04 
0 . 312015700421D-03 


ERROR 

-0. 17685419 1014D- 12 
0 . 496907202885D-11 
-0.713074044256D-11 
-0. 866283433876D-11 
0.698968799506D-10 
-0.9 395520472 74D-10 
-0 . 364804006447D-09 
0.3 87429 1026 88D-08 


ERROR 

-0 . 128776609562D-01 
0 . 3Q2541807675D-01 
-0. 55832 765 399 ID-01 
0.149063705060 
-0.319785257172 
0.348765545461 
0.439916989730 
-4.27885454209 
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An evaluation of the previous tests must take into account these 
three areas of performance: 1) accuracy of solution, 2) cost and 

3) ease of programming* To make a judgment in the first category we 
must consider Sections 3 and 5* Only in these sections were the true 
solutions known. Table 3*1 shows that for the type of matrix in question 
neither method of solution has a decided advantage over the other. Both 
have almost the same limiting values for noise levels and condition 
numbers. For the matrix in Section 5, it was found that Gaussian 
elimination worked better than the optimization method in every case. 

It was found that for identical problems the optimization method 
cost more than three times the other method. Not only is the cost 
prohibitive, but as was mentioned earlier, DFMFP had a tendency to go 
into underflow or overflow. This could be avoided if a multiplying 
factor was introduced in the subroutine FUNCT. However, it was not 
clear what this factor should be and thus necessitated trial-and-error . 

Of course, this is not only inconvenient to the programmer but also 
adds to the already high cost of the optimization method, It can now 
be said that in the solution of a matrix equation, this optimization 
technique holds no advantage over the usual Gaussian elimination and, 
in fact, is inferior in all categories considered. 
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